


X) 



Oh' 

o 



X 



Submitted to the Astrophysical Journal 

Dynamics of Line-Driven Winds from Disks in Cataclysmic Variables 
II. Mass Loss Rates and Velocity Laws 

Achim Feldmeier 
Department of Physics & Astronomy, University of Kentucky, Lexington, KY 40506-0055, USA, E-maih 

achimOpa . uky . edu 

and 



Isaac Shlosman 
Department of Physics & Astronomy, University of Kentucky, Lexington, KY 40506-0055, USA, E-mail: 
^ ' shlosmaji@pa.uky.edu 

,_- . and 

Os 

Peter Vitello 

^ ■ Lawrence Livermore National Laboratory, L-14, P. O. Box 808, Livermore, CA 94550, USA, E-mail: 

O^ , vitello@spaniel.llnl.gov 

(NJ ■ ABSTRACT 

o 
(^ 



We analyze the dynamics of 2-D stationary, line-driven winds from accretion disks in cata- 
clysmic variable (CV) stars, by generalizing the Castor, Abbott and Klein (CAK) theory. In the 
first paper, we have solved the wind Euler equation, derived its two eigenvalues, and addressed 
the solution topology and wind geometry. Here, we focus on mass loss rates and velocity laws of 
the wind. 

^ [ We find that disk winds, even in luminous novalike variables, have low optical depth, even in 

J^ ' the strongest driving resonance lines. This suggests that thick-to-thin transitions in these lines 

occur in the wind. For disks with a realistic radial temperature law, the mass loss is dominated by 
gas emanating from the inner decade in radius. The total mass loss rate associated with the wind 



^ • from a disk of luminosity 10 L© is ~ 10~^^ M0 yr"^, or 10^^ of the mass accretion rate. This is 

^ one order of magnitude below the lower limit obtained from fitting P Cygni line profiles using 

kinematical wind models, when the ionizing flux shortwards of the Lyman edge is supressed. The 
difficulties associated with such small mass loss rates for line-driven winds from disks in CVs are 
principal, and confirm our previous work on this subject. We conjecture that this issue may be 
resolved by detailed non-LTE calculations of the CAK line force within the context of CV disk 
winds, and/or better accounting for the disk energy distribution and wind ionization structure. 
We find that the wind velocity profile is well approximated by the empirical law used in 
kinematical modeling. The acceleration length scale is given by the footpoint radius of the wind 
streamline in the disk. This suggests an upper limit of ~ 10 r^d to the acceleration scale, which 
is smaller by factors of a few as compared to values derived from line fitting. 

Subject headings: accretion disks — cataclysmic variables — hydrodynamics — stars: mass-loss 
— stars: winds 



1. Introduction 

Line-driven winds (hereafter LDWs) are expected around luminous objects whose spectra peak in the 
UV, such as OB stars and accretion disks, stellar and galactic. Feldmeier & Shlosman (1999, Paper I) have 
investigated a 2-D analytical model of LDWs from disks in cataclysmic variables (CVs) characterized by a 
large mass transfer rate from the secondary to the white dwarf. Such CVs, i.e., novalike variables and dwarf 
novae in outburst, show clear signs of outflows driven by radiation pressure (Paper I and refs. therein). 

Recent numerical simulations of time-dependent 2-D disk winds by Proga, Stone, & Drew (1998; PSD 
hereafter) largely confirmed the previous kinematical studies by Shlosman & Vitello (1993) and Knigge, 
Woods, & Drew (1995), and delineated a number of empirical relationships which require further physical 
explanation. Paper I has addressed the 2-D geometry of the wind streamlines and the topology of solutions 
to the wind momentum equation. In particular, a comparison with the PSD model was made as well as a 
comparison with the 1-D LDWs from OB stars. 

The main results of Paper I are as follows. Firstly, the solutions to the wind momentum equation are 
characterized by two eigenvalues, the mass loss rate and the flow tilt angle with the disk, in the presence 
of a realistic radiation field above the disk. The existence of the second eigenvalue is a reflection of the 
multi-dimensional nature of a disk wind. The wind itself appears to be collimated to a certain degree, i.e., 
the wind coUimation angle with the rotation axis (semi-opening angle) is about 10° for the wind launched 
within 4 white dwarf radii, and about 30° for the outer wind, for the Shakura-Sunyaev (1973; SHS hereafter) 
disk. Furthermore, the wind coUimation depends solely on the radial temperature stratification in the disk, 
unless there is an additional degree of freedom such as a central luminosity associated with nuclear burning 
on the surface of the white dwarf. The above degree of coUimation for disk winds in CVs should be taken 
with caution at large distances from the disk. 

Secondly, a major distinction between stellar and disk LDWs is the existence of maxima in both the 
gravity and the radiation flux above the disk. This behavior of gravity and radiation flux results in profound 
topological differences in the solutions to the stellar and disk wind momentum equations. We find that two 
distinct regions of disk wind exist, the inner and outer winds. The critical point of the outer wind lies close 
to the disk photosphere and to the sonic point, upstream of the top of the gravity 'hill'. This proximity of 
the critical and sonic points is typical of LDWs from O stars as well. On the other hand, for the inner wind, 
the critical point lies far away from the sonic point, beyond the gravity hill. 

Observationally, the mass loss rates from CV disk winds are poorly constrained. This is mainly a 
consequence of uncertainties in the ionizing fluxes from different components in the system, and, therefore, 
of the ionization stratification of the wind. Neglecting the boundary layer and assuming a local blackbody 
radiation from the disk, Vitello & Shlosman (1993) and Knigge et al. (1995) find wind mass loss rates 
^ 1% of the accretion rate, by fltting observed P Cygni line proflles. For a system of luminosity 10 Lq, this 
corresponds to a mass loss rate of about 2 x 10^^" Mq yr~^. 

These mass loss rates are upper limits for the following reason. The continuum radiation field shortwards 
of the Lyman edge is often found (Polidan, Mauche, & Wade 1990; Long et al. 1991, 1994; van Teeseling, 
Verbunt, & Heisc 1993; Knigge et al. 1997) to be highly supressed compared to blackbody emission or stellar 
photospheric fluxes, or is even absent. This drastic reduction in ionizing flux allows to reduce the electron 
density, and, therefore, M, while maintaining the same degree of ionization in the wind. 

On the other hand, a reasonable lower limit to the wind mass loss rate from luminous CV disks was 
found by Prinja & Rosen (1995). They argued that the product Mq, where q is the ionization fraction of 



Mqjt^'^ for ten dwarf novae and novalike variables with high- 
resolution lUE spectra (see also Mauche & Raymond 1987; Hoare & Drew 1993). This results in a lower 
limit of M ^ 5 X 10~^'^ Mq yr~^. Note that this value is still model dependent to some degree, since Prinja 
& Rosen assumed a constant ionization fraction throughout the wind, and a purely radial flow. 

In this paper, we focus on the mass loss rate and the velocity law of disk LDWs in CVs. We employ 
the main assumption of Paper I that the wind streamlines are contained in straight cones, whose coUimation 
angles are eigenvalues of the flow. This approximation is acceptable close to the disk (e.g., as shown by 
PSD) and allows us to make a meaningful determination of the mass loss rate (which is constrained close 
to the disk photosphere as well) , and of the initial velocity law. A similar conclusion was reached from line 
fitting using 2-D wind kinematics (Vitello & Shlosman 1993; Knigge et al. 1995). On the other hand, we 
expect both the centrifugal forces and the polar component of the line force to bend the streamlines at large 
distances from the disk, and to influence the terminal velocity of the flow. 

This paper is organized as follows. Section 2 reviews the mass loss rates derived from CAK theory. 
Section 3 shows that CV disk winds are expected to have low optical depth even in strong lines. Sections 4 
and 5 derive the mass loss rates and velocity laws for the disk wind model. Sections 6 and 7 discuss and 
summarize our results. 



2. CAK mass loss rates for stars and CV disks 

The CAK line force for stellar winds is fully determined by two parameters, the power exponent a 
and the mass absorption coefficient of the strongest driving line, kq (Paper I; Puis, Springmann, & Lcnnon 
1999). Instead of using kq directly, CAK parameterize the line force per unit mass as gi, = M(t)ge, where 
M{t) = kt'" is the so-called force multiplier, and ^c is the force due to electron scattering. The optical 
depth t refers to a line with k = Ce, where a^ is the electron scattering coefficient. The parameter k is given 
in terms of kq by 

)b = iM^f^V"", (1) 

1- a c \aej 
where T{a) is the complete Gamma function, and Wth is the thermal speed of carbon ions (CAK) or hydrogen 
(Abbott 1982). The above parameterization for M{t) as a power law, however, does not account correctly 
for optically thin winds, where the force multiplier saturates at some value Mj^a^it) = Q ~ 2, 000 for a gas of 
solar composition (Abbott 1982). This particular value of Q was determined for O stars. However, already for 
winds from B stars near-the-main sequence, modeling of X-ray spectra suggests that mass loss rates inferred 
from the standard theory can be substantially underestimated (Cassinelli 1994; Cassinelli et al. 1994). In 
terms of Q, kq is given by (Gayley 1995), 

^^ = r(,)-T^ Q. (2) 

Inserting the eigenvalue E^ from eq. (1.8b; hereafter tables, figures and equations from Paper I have a roman 
prefix T) into (1.6), the CAK mass loss rate from O star winds takes the compact form. 



1-r. 



M-T^(^) 5- (3) 



This expression is valid as long as QT/{1 — T) > 1. Otherwise, the gravity prevails and no wind solution 
exists. For the above value of Q, and assuming a = 2/3 and an Eddington factor F = 0.5 of an O supergiants. 



eq. (H) gives M ^ 90 L/c^. The mass loss rate from CAK theory agrees thus well with the estimate from 
the single scattering limit, M = (c/voo) {L/c^) ~ 100 L/c^, for typical terminal speeds v^c ^ 3,000 kms~^. 
However, agreement between both mass loss rates is solely due to the fact that V is close to unity. 

If, alternatively, F <C 1, the CAK mass loss rate falls well below the single scattering limit. As we show 
in the next section, the mass loss of a disk LDW is given again by eq. (||), up to correction factors of order 
unity. Even for the brightest CVs, i.e., novalike variables and dwarf novae in eruption, which experience 
LDWs, QT ^ 1 (applying the O star value of Q). Hence, M ~ L/c^ from eq. (^, whereas the single 
scattering limit gives M ^ 60 L/(?, for v^o ~ 5, 000 kms~^. Since in thin LDWs the probability for a photon 
to be scattered by a line is less than unity, the estimate from the single scattering is way too high (see also 
Puis, Springmann, & Owocki 1998). 



3. Disk wind optical depths 

A number of fundamental differences exist between stellar and disk LDWs in CVs, some of them dis- 
cussed in Paper I. Here we show that optical depths for CV winds are more typical of thin winds, e.g., of 
B stars near-the-main sequence, than of more extensive supergiant winds. 



3.1. Disk virind geometry and radiation field 

The disk wind geometry is described in Paper I, and we repeat here only the essential assumptions and 
make necessary definitions. A flow streamline is a helix which is contained within a straight cone (Fig. 1). 
The footpoint radius of a streamline in the disk is ro, the tilt of the cone with the disk is A, and x is the 
distance along the cone. We neglect pressure forces, and assume that the azimuthal velocity is determined 
by angular momentum conservation above the disk, and by Keplerian rotation within the disk. The only 
remaining velocity component is v^, which points upwards along the cone. We introduce a normalized 
coordinate X = x/rp. The velocity V is normalized to the local escape speed, and the flow acceleration 
becomes W' — 2V dV/dX. Lastly, we introduce the radiation flux F above the disk and the flux F normalized 
to the footpoint flux of the streamline, as well as their projected counterparts F^ and F^ along the wind 
cone (Section 3.2 of Paper I). 



3.2. Semi-transparent Vifinds from disks in CVs? 

In this section, we make use of a simplifled Euler equation for the disk LDW to show explicitly that low 
F factors in CV disks imply low optical depths in the wind. By doing so, we neglect factors of order unity 
from the angle integration in eq. (1.2), in replacing r-y by t^. For simplicity, only the disk with F (x r^ 
(hereafter 'Newtonian') is considered. The line force is then 

cTeFo,{ro,X) ge-Fz(r-o,0) ^ 

Sl = M[t^) = Fx(ro,X)M{tx), (4) 

c c 

The Euler equation for the disk wind, in the limit of zero sound speed, and neglecting the force due to 
electron scattering because of small F above the disk, becomes (with M^j being the white dwarf mass). 



, G Afwd 

«^^x = -2— 

'o 



g{X)-T{ro)F4ro,X)M{t^) , (5) 



where the effective gravity g was defined in (1.9), and wc introduced 

. ^4?!^. (6) 

For the Newtonian disk, F = cTg id/[47rcG'Mwdln(rd/r„d)] becomes independent of tq. Here, L^ is the 
disk luminosity and r^ and r^d are the outer and inner disk radii, r„d being the white dwarf radius. Using 
typical parameters for novalike CVs, L^ = 10 Fq, M^d — 0.7 Mq, and rd/r^d — 30, one has that F ~ 10^"*. 

For a stationary wind solution to exist, the right-hand-side of eq. (g|) has to be positive. This poses a 
constraint on M{t) and therefore on t. Namely, the maximum of reduced gravity g lies between 2/(3V3) — 
0.38 for A = 90° and 4/27 ~ 0.15 for A = 0°. Since F^ is of order unity (see Fig. F5), M{t) approaches its 
maximum value, Q, in regions of large gravity, and stays constant thereafter. In other words, because F is 
so small for CV disks, the wind solution barely 'makes it' over the gravity hill. 

This saturation effect in M{t) happens when the strongest driving line in the wind becomes optically 
thin, at about t ^ 10^^. If this thick-to-thin transition occurs before or at the critical point of the flow, 
the wind solution is lost. This possiblity cannot be excluded in our model due to the rapid change of the 
velocity gradient (and hence of t) in the vicinity of the critical point. The consequences of this effect on the 
feasibility of FDWs are discussed in Section g. 

The situation is fundamentally different for dense winds of O stars, where the Euler equation reads 

vv' = -^^{l-T[M(t) + l]). (7) 

Assuming F > 0.1 for O stars, there is a wide range in M{t) for a stationary solution to exist, namely from 
~ 10 — 2,000. The highest mass loss rate solution (hence, the slowest wind) is characterized by the lowest 
allowable M{t) ^ 10. The permitted range in M{t) corresponds to an even wider range in t (because a < 1), 
t ^ 10^^ — 10"'^. We further quantify these arguments in Appendix A. 

Our estimate for t in CV winds contradicts the claim by Murray & Chang (1996) that the optical depth 
parameter t is similar for CV disk and O star winds. 



4. Disk wind mass loss rates 

4.1. Vertical disk wind 

The mass loss rate of the vertical wind above an isothermal disk is determined by the eigenvalue Ec in 
eq. (1. 18b). From eq. (1.17) one finds, 

l-Q 

M^^(^-^h^ " d(^)\ (8) 

1-a \^87rGAfwdy \c^ J ^ ' 

The dimensionless constant D is given by 

i. = /^(^--^ff'°y, (9) 

ri 

where ri and r2 are the inner and outer radii of the wind base in the disk, respectively. Making the plausible 
assumptions that r^ 3> r^d, ^2 2> ^i, and r2 ~ r^ (both of the latter radii are determined essentially by 



the temperature dropping below 10"* K), we estimate D ~ 1 for a = 1/2 and 2/3. Introducing the disk 
Eddington factor 



eq. (H) becomes 



(10) 



M-^Y^\^Q^^] TT- (11) 




Up to the correction factors due to different gravity (and geometry), this equation is identical to the CAK 
mass loss rate from a point star, eq. (y). Note that a disk wind is more efficient in carrying mass loss than 
an O star wind by a factor (3\/3/2)*-^~"-*/", which is due to the lower gravity potential well. Re-writing 
eq. (O) as M = N L^/c^ (where the coefficient N depends on Ld), and using the relevant parameters for 
novalike CVs introduced above, gives TV ~ 2 and a disk mass loss rate M ~ 10""'^^ M© yr~^, for (5 = 2, 000. 



4.2. Tilted Disk Winds 

A more realistic picture of disk mass loss consists of a tilted wind from a disk with radial temperature 



stratification. From eq. (1.17), using eq. (10), 



where 



M = ^(^Qr,) " i.(..)^, (12) 



" -m. PF,(r,0)]'/"/ r ^pF,(r,0)]'/", (13) 



and TTir is normalized in units of mass loss from a vertical wind; (m) is then a normalized mass loss 
rate per dr, averaged over the wind base. Using values of rhr from Table I.l, we estimate (to) ~ 1.2 
for a = 2/3, and (to) = 2.2 for a — 1/2, both for SHS and Newtonian disks. For the latter, one has 
D = ln(r2/ri)/[ln(rd/r„d)]^^"- This expression is roughly correct also for the SHS disk, where D cannot 



be calculated analytically. Inserting the values of (to) and D into (12) one finds, for a = 1/2 and 2/3, and 
for SHS and the Newtonian disks, assuming once again typical parameters for novalike CVs and Q — 2, 000, 
that M ~ 10^12 Mgyr-i. 

Interestingly, the mass loss rates for a vertical wind above an isothermal disk and for tilted winds above 
disks with temperature stratification are very similar. This means that the disk mass loss rate is only a weak 
function of the tilt angle, as long as the latter is aligned with the radiation flux. For the disk types used in 
the present work, this range encompasses A ~ 50 — 90°, according to Table I.l. 

It is readily shown that the mass loss rate for a LDW due to a single, optically thick line is roughly 
Ld/(? ■ The above N = 2, therefore, implies that only few lines become optically thick in the present CV 
disk wind model. Unlike for disk winds, in O star winds of the order of 100 lines become optically thick 
according to eq. (||). 

Note that in Paper I, eigenvalues E were derived without including the saturation of the force multiplier 
at Minax(0 = Q (i-6-, without applying the exponential line-list cutoff of Owocki, Castor, & Rybicki 
The above mass loss rates are, therefore, upper limits. 



Finally, we derive the dependence of the mass loss rate dM on tq. According to (1.17), this relation is 
determined by the disk temperature stratification and the run of Ec with rg. For an isothermal disk, E^ 
was found to be independent of rg, hence dM / dr^ oc r^ "-!/"_ Such an unrealistic growth of dM /dr^ with 
radius is a consequence of the increase of radiation energy with rg . Alternatively, for the Newtonian disk, 
fitting Ec{ro) from Table I.l with a power law, one finds approximately dM/drg oc r^ ' for rg ^ 5r„d- 
The total disk mass loss rate scales therefore roughly as M c>c lnr2. For the SHS disk, approximately 
dM(ro)/dro ex Tq ' for rg ^ 5rwd- At smaller rg, the radial dependency is weak for both types of disks. 
Therefore, the mass loss rate M is centrally concentrated, and more so for SHS disks. 



5. Velocity laws for disk winds 

We discuss the wind velocity law by solving the Euler equation, first neglecting, then accounting for 
ionization effects in the fiow. Only the flow above the SHS disk is analyzed in this section. 



5.1. Solutions to the algebraic Euler equation 

The geometrical expansion terms, discussed in Section 5.1 of Paper I, introduce a velocity dependency 
into the wind Euler equation. However, we find that these terms leave the mass loss rates in the model 
practically unchanged, and increase the wind velocity by at most 10%. The reason for this is the small angle 
dispersion in the wind streamlines. We, therefore, omit geometrical expansion terms here. 

Without explicit dependence on velocity, the Euler equation becomes purely algebraic. Figure 1.8 
displays the solutions W'{X) of the Euler equation (1.20) for a tilted disk wind, at different rg. The velocity 
field above the disk is obtained by integrating W' , and is displayed for a number of streamlines in Fig. g, for 
the SHS disk and a = 2/3. As shown in Paper I, the LDW velocity law is not a monotonic function of X. 
Clearly, the deceleration regime (marked with a cross in Fig. 0) is unimportant for streamlines starting at 
large rg, because it lies at large X where the fiow moves much faster than the local escape speed. The inner 
streamlines do show a more pronounced kink. 

The wind terminal velocities are found to be independent of rg, for a — 1/2, and linearly dependent 
on rg for a = 2/3, growing by a factor of five over the whole disk. However, the above examples may 
be of academic interest only, since the wind is expected to go through an optically thick-to-thin transition 
beforehand. Hence the actual observable terminal velocities may be smaller. Note that we presume that the 
wind geometry stays unchanged at all X, and neglect the saturation effect in the force multiplier at very 
small optical depths in the wind. Both effects can influence the asymptotic wind dynamics. 

The velocity profiles in Fig. are well approximated by the empirical velocity law used by Shlosman & 
Vitello (1993) in fitting the line profiles of novalike CVs, 

.-. (^/^--)' (14) 

""-""^i + ix/x^^^r- ^^^> 

Here, ATacc is the acceleration length scale along the wind cone. We find the best fits to Fig. for /? = 1.5 — 1.9. 
Vitello & Shlosman quote rather similar values of 1.3 — 1.5. Furthermore, ATacc ~ 1 in Fig. g, which means 
that the footpoint radius rg sets the acceleration length. The reason for this is that Xacc is determined by 
the effective gravity and disk radiation field, i.e., by the auxiliary functions g and / (Paper I), which change 



on length scales ~ tq. For this reason, we do not expect Xacc to depend on a. The Xacc from observed 
P Cygni line profiles are found to be somewhat larger, namely in the range 1 — 10 (Hoare & Drew 1993; 
Vitello & Shlosman 1993), depending on individual objects. 



5.2. Effects of ionization stratification in the wind 

The major concern of the present model is the small optical depth in the wind, where the latter appears 
to be semi-transparent even in the strongest driving lines, when Q — 2, 000 is used. A possible resolution 
of this problem may be related to our neglect of ionization structure in the wind, which is expected to 
lead to a shallower velocity law and smaller velocity gradients. To proceed, we parameterize the ionization 
stratification in the simplest possible way, by introducing the 5{> 0) parameter customarily used in stellar 
wind theory (e.g., Abbott 1982), namely 

5L « t-"C\ (15) 

where S, = J/f^o is the ionization parameter, with J being the frequency-integrated mean intensity. Higher 
ionization stages typically harbor fewer lines than lower ionization stages, and therefore lead to a smaller 
line force. Typical values for O stars are 5 <^ 0.1 (Abbott 1982; Puis 1987; Pauldrach et al. 1994), but 
values as large as 5 — 0.7 have been suggested recently for winds at low effective temperatures of 8,000 K 
(Kudritzki et al. 1998). 

To include the 5 correction term in the present model, we calculate first dJ from a disk annulus of width 

\/ q^ + r"^ + z"^ — 2rq / Arq 

(g2 + r2 -f z2)2 - {2rqf ^ \J^ 

where / is the isotropic intensity, and E2 is the complete elliptical integral of the second kind (Abramowitz 
& Stegun 1965). For the temperature stratifications of interest, this expression for dJ cannot be integrated 
analytically over q to give J. Resorting to a 1-D numerical integration. Fig. shows J and F as function of X 
for the SHS disk at some representative ro. Importantly, for the outer wind, the critical point lies upstream 
of the maximum of both J and F. Since J increases along the streamline while the electron density drops, ^ 
increases all the way to the critical point, as in O star winds. The line force at the critical point is, therefore, 
smaller for S > than for ^ = 0, and the same is true for the mass loss rate. Alternatively, in the unlikely 
situation that the wind recombines along the streamline, the driving force as well as the mass loss would 
increase, because of a larger number of metal transitions. 

In order to understand the effect of 6 on the wind velocity law, we include the (5-correction (or ionization 
stratification in the wind) above the critical point but neglect it below the latter, thereby leaving the mass 
loss rate unaltered. Since the J-term introduces a dependence of the line force on W, the solution to the 
Euler equation has to be iterated until convergence is achieved. 

We find that, assuming a = 2/3 and S = 0.2, terminal speeds decrease by a factor of ~ 2 for outer 
regions of the SHS disk, whereas the optical depth t increases by a factor of ^ 4, thereby pushing the 
solution further away from the cutoff at i ^^ 10~^. For inner disk regions, the effect of S on v^a and t is less 
pronounced. The somewhat ambivalent conclusion is, therefore, that 5-terms may, as desired, raise t by a 
factor of a few, but, at the same time also can lead to an unwanted reduction in terminal speeds. 



dJ{r,z) ^ 4:I{q,0)qdqz ^.^ , ^2 , — 2\2 a7Z:T2 -^2 ( -77 ^^2 , ,2 ) ' i^^) 



6. Discussion 

In this section, we compare the mass loss rates from our model with values quoted in the literature, 
which are typically estimated from P Cygni line fits or from dynamical wind modeling. We also mention 
briefly some processes neglected in this work, which may have an effect on M . 



6.1. Comparison with M from kinematical wind modeling 

The most reliable estimates for mass loss rates from cataclysmic variables are so far from P Cygni line 
fits. For novalike variables with parameters similar to those considered here, Vitello & Shlosman (1993) 
and Knigge et al. (1995) find a lower limit of ~ 10~^ to the ratio of mass loss-to-accretion rates. Below 
this value, the observed line profiles cannot be reproduced because the wind becomes over-ionized. This is 
about two orders of magnitude larger than M derived in Section ^, when Q = 2,000 is used as in O stars. 
It is important, however, that in the line fitting a highly idealized blackbody disk spectrum was used. 
This is a clear overestimate of the ionizing radiation flux shortward of the Lyman edge (Polidan et al. 1990; 
Long et al. 1991, 1994; van Teeseling et al. 1993; Knigge et al. 1997). For a more realistic EUV radiation field 
which accounts for the Lyman continuum cutoff, the wind mass loss rate can be reduced, while maintaining 
the same ionization parameter. We have been able to produce acceptable P Cygni line profiles down to 
M ~ 10""'^^ Mq yr~^, for the present system parameters. Essentially all the carbon was then in the form of 
C IV, due to Auger ionization of C ii by X-rays. 



6.2. Comparison with M from dynamical wind modeling 

Next, we compare the mass loss rates from our analytical wind model with those from dynamical 
simulations of disk winds. The only realistic dynamical modeling of CV disk winds attempted so far was 
performed by PSD. We have provided a general comparison between our models in Paper I, and turn here 
to mass loss rates. 

Using k = 0.2 and a = 0.6, PSD find for the SHS disk mass loss rates of M = 5 x W'^'^MQyT''^ and 
M = 5 X lO^^^Moyr"^, corresponding to disk luminosities oi L = 8 L0 and i = 24 L0, respectively. It 
is clear that such an increase by a factor of 100 in M cannot be understood from the simple CAK scaling, 
where M ex L^^". Instead, this strong dependence on L^ is a consequence of the optically thick-to-thin 
transition in the disk LDW when QT ^ 1. The reason for this is that PSD apply the exponential line-list 
cutoff (Owocki et al. 1988), due to which the force multiplier M{t) reaches a maximum value in optically 
thin flow regions. This suppresses the mass loss rate as compared to the case of a pure power law force 
multiplier. 

Even, for a low-luminosity disk, PSD find that large M can be driven when a is as large as 0.8. However, 
PSD assume the same value of fc = 0.2 for all a. As is evident from Fig. 4 of Gayley (1995), k drops by a 
factor of 10 when a increases from 0.6 to 0.8, leading to a very similar M in both cases. This is a consequence 
of Q r ~ 1 for a disk wind in eq. (p2). Contrary to this, for O star winds, QT ^ 1, hence M depends strongly 
on a. 

For a more luminous disk with L^ = 24 L©, we estimate M ~ 4 x 10"^'^ Mq yr^^ from eq. ( [l2|) when a 
value of D is used characteristic of the narrow wind base of PSD. The agreement with Af = 5 x 10~^^ M© yr~^ 
as found by PSD is very good. 
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6.3. The mass loss paradigm in disk winds from CVs 

Taken at face value, the modified CAK theory of LDWs from CV disks predicts surprisingly low mass loss 
rates when a line force parameter Q = 2, 000 is used as in O stars. The calculated rates of ^ 10^^^ Mq yr^^ 
for Ld = 10 Lq (or even lower, when the saturation of the force multiplier is accounted for), imply that LDWs 
will have a thick-to-thin transition in the strongest driving lines, and will have difficulties in reproducing the 
observed line profiles. What are the possible solutions to this problem? 

Firstly, higher values of Q are the most obvious way to increase M in the wind. The value of Q is 
highly uncertain. Especially, it is unclear how the different spectral shape of disk radiation and its effect 
on the non-LTE electron occupation numbers in ions will modify Q as compared to its O star value. Such 
calculations have never been attempted and are clearly beyond the scope of this paper. Theoretically, the 
value of Q can be significantly larger than for the O stars. 

Analysis of O star LDWs revealed Q Gamma to be the key parameter determining the mass loss rate 
in the wind (Gayley 1995). Namely, for QF < 1, the wind ceases to exist. However, for O stars QF > 100, 
producing a kind of a 'safety belt'. The situation is different for winds from luminous CV disks, where 
QT ~ 1, if the O star value for Q is used (this work and PSD). This was already foreseen by Vitello & 
Shlosman (1988), which claimed that the disk winds appear to be much more restrictive than O star winds. 

Secondly, the CAK line force does not account for the ionization stratification in the wind. So a situation 
can arise where the wind is driven by the part of the disk which is UV bright, while the ionization is controlled 
by the central disk region, e.g., the boundary layer, which is X-ray bright. This will be especially pronounced 
in the presence of a Lyman continuum cutoff in the UV source. Under these circumstances, due to the wind 
ovcr-ionization, the line force will be reduced below its CAK value (MacGregor, Hartmann, & Raymond 1979; 
Fransson & Fabian 1980) and hence the velocity gradient as well. A shallow velocity law will increase the 
overall optical depth in the wind and push the thick-to-thin transition downstream, away from the critical 
point. It is important in this context that the minimal changes in the P Cygni line profiles in eclipsing 
novalike CVs require very shallow velocity profiles in the wind (e.g., Shlosman, Vitello, & Mauche 1996). 

Additional processes can affect the mass loss rates, however to a lesser degree. Firstly, the assumption 
that each wind parcel is confined to a straight cone is a major limitation of the present model. For bended 
streamlines, a better alignment with the radiative flux vector seems possible, which should increase the mass 
loss. Secondly, Pauldrach & Puis (1990) find a sudden increase in stellar mass loss rates of B supergiants 
when the Lyman continuum becomes optically thick at the wind base, e.g., when the effective temperature 
drops below a certain threshold or if F reaches some critical value. This induces a shift to lower wind 
ionization, which in turn increases the mass loss (more driving lines), increases the Lyman jump even more 
and so on, until a stable situation with high M and low wind ionization is reached. Thirdly, Q depends 
linearly on the wind mctallicity, and enhanced metal abundances could lead to a larger QT, which governs 
the mass loss. This is significant in particular to the He abundance due to the efficient convection in the 
low-mass secondary stars, but may be relevant to metals as well. 



7. Summary 

The focus of this second paper on line-driven winds from accretion disks in cataclysmic variables is on 
the theoretical estimates of the mass loss rates and wind velocity laws. Our results arc as follows. 

The derived mass loss rate using the modified CAK formalism appears to be substantially smaller than 
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that inferred from P Cygni line fits, even with supressed Lyman continuum, and more so when the saturation 
effect in the line force multiplier is included. Yet this rate is in a good agreement with results from time- 
dependent, 2-D dynamical wind simulations by PSD. The reason for low mass loss rates is that the key 
parameter controlling LDWs, QT ^ 100 for O star winds, is only ^ 1 for disk winds, when O star values for 
Q are used. Some potential resolutions to this problem were proposed. 

We find that the mass loss is dominated by the inner decade in disk radii. For the Shakura & Sunyaev 
and for Newtonian disks, the mass loss per unit radius is roughly uniform out to five white dwarf radii, and 
drops oc r~^'^ and r~'^'^ at larger r, respectively. 

Due to their low mass loss rates, CV winds should experience a thick-to-thin transition even in the 
strongest lines. CV winds should, therefore, resemble more closely winds of B stars near-the-main sequence 
than of O supergiants. 

The wind velocity profiles show a slowly accelerating flow, with a characteristic acceleration length given 
by the footpoint radius of the streamline in the disk. Fitting the observed line profiles using kinematical 
models suggest even slower accelerating winds. The observable terminal velocity of the wind is associated 
with the thick-to-thin transition in the driving lines. Due to this latter fact, and due to uncertainties in the 
ionization stratification and in the anticipated streamline bending at large radii, the actual wind terminal 
velocity is poorly constrained in our model. 

We are grateful to Jon Bjorkman, Rolf Kudritzki, Chris Mauche, Norman Murray, Stan Owocki and 
Joachim Puis for numerous discussions on various aspects of line-driven winds. I.S. acknowledges hospitality 
of the IGPP/LLNL and its Director Charles Alcock, where this work was initiated. This work was supported 
in part by NASA grants NAG5-3841 and WKU-522762-98-06, and HST AR-07982.01-96A (I.S.), and was 
performed under the auspices of the US Department of Energy by LLNL user contract number W-7405- 
ENG-48 (P.V.). 



A. Disk Eddington factors required for disk line-driven winds 

Using the eigenvalues E^ from Paper I, one can further quantify the disk Eddington factors required to 
launch a line-driven wind which is optically thick at least up to the critical point. Consider first a vertical 
wind above an infinite, isothermal disk. The Euler equation (g) becomes, for F^ = 1 (the flux component 
along X, normalized to the footpoint flux), 

W'^-g + Tiro)M(U). (Al) 

From eq. (1.18) at the critical point, 

r(ro) - , ^^^^,, . (A2) 

^°' (l-a)Me(t) ^ ' 

Since the maximum of M{t) for a gas of solar composition is ^ 2, 000 (Abbott 1982), one has that r(ro) > 
4 X 10^"* for all vq. Next, for the more realistic case of a tilted wind above a disk with radial temperature 
stratification, the simplified Euler equation (||) reads 

W ^-g + r{ro)F,M{t,). (A3) 



12 



Using W^ = {a/1 — a) gc and assuming gc ~ Xc for critical points close to the disk, 

r(ro) = ^ . (A4) 

{l-a)F^{ro,X,)M,{t) 

According to Fig. 1.5, F < 4, hence r(ro) ^ 10^'* is required for the wind to be optically thick at the critical 
point. 
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Fig. 1. — Adopted flow geometry for a CV disk wind. The streamlines are helical lines, and are assumed to 
lie on straight cones. 
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Fig. 2. — Velocity profiles for the SHS disk with r^ — 30rwd, for a — 2/3, and at footpoint radii 3, 7, 10, 
ISr^d (from bottom to top). The dotted lines are fits using the velocity law eq. ([l4|), with the values of 
(Xacc;/3) given in brackets. The velocity V is in units of the local escape speed at vq. For the tilt angles, 
eigenvalues Ac from Paper I were used. Triangles mark the critical point, crosses the point where the wind 
starts to decelerate. 
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Fig. 3. - Normalized, projected flux (full line) and mean intensity (dashed line) above the SHS disk 
(rd = 30rwd), along a ray with footpoint tq = lOrwd and Ac — 58°. 



